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^ SUMMARY 
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source  files  and  their  use  on  the  A RL  computer,  and  suggestions  on  how  to  implement  the 
language  on  other  computers. 
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1.  INTRODUCTION 


Since  documenting  the  simulation  language  CSMP-IO(ARL)  (see  Ref.  1 and  2),  there  have 
been  a number  of  developments  necessitating  the  following  postscript.  These  developments, 
described  below,  include  (a)  the  provision  of  an  alternative  integration  method,  (b)  compilation 
using  the  DEC  FORTRAN- 10  compiler,  and  (c)  a number  of  improvements,  especially  to  the 
output  program.  Some  of  these  developments  are  included  in  Reference  3,  but  are  only  described 
very  briefly.  Developments  (a)  and  (b)  have  resulted  in  four  versions  of  the  modelling  program 
BOMMP.  In  order  to  avoid  unnecessary  repetition  of  stored  information,  the  source  files  have 
therefore  been  suitably  reorganized.  The  location  of  the  necessary  source  files  at  ARL  and 
details  on  how  to  obtain  running  versions  on  the  ARL  computer  are  given  here.  Following 
interest  shown  in  using  the  language  on  other  computers,  the  uses  made  of  a small  number  of 
subroutines  written  in  the  DEC  machine  language  MACRO-10  have  been  investigated.  Descrip- 
tions of  the  operation  performed  by  each  of  these  routines  are  therefore  provided  here  together 
with  suggestions  on  action  to  take  when  alternative  subroutines  are  unavailable. 

The  aircraft  arresting  gear  problem  of  Reference  4 is  used  to  fully  demonstrate  the  new 
operating  procedures  and  improvements  to  the  simulation  language  (see  Appendix  A). 

In  References  1 and  2,  when  showing  terminal  messages,  contrasting  type  faces  were  used 
to  distinguish  between  messages  typed  by  the  user  and  those  typed  by  the  appropriate  computer 
program.  Once  familiar  with  CSMP-I0(ARL),  the  above  style  of  presentation  becomes  un- 
necessary. Hence,  to  avcid  high  typesetting  costs  for  this  publication,  the  actual  terminal  printout 
is  used  directly  in  Appendix  A. 

2.  THE  MODELLING  PROGRAM  ‘BOMMP’ 

2.1  Improvements 

The  improvements  described  below  in  Sections  2.1.1  and  2.1.2  have  resulted  in  four  versions 
of  BOMMP  which  are  summarized  in  Table  1.  It  should  be  noted  that  Versions  2 and  4 may 
also  be  compiled  using  the  FORTRAN  IV  (i.e.  F40)  compiler  as  well  as  the  FORTRAN-10 
(i.e.  F10)  compiler.  Whichever  compiler  is  used  however,  Versions  2 and  4 do  not  provide 
expandable  arrays. 


TABLE  1 

Versions  of  BOMMP 


Version 

Compiler 

Adams 

Integrator 

Expandable 

Arrays 

1 

F40 

no 

yes 

2 

FIO 

no 

no 

3 

F40 

yes 

yes 

4 

F10 

yes 

no 

2.1.1  Alternative  Integration  Method 

An  Adams  predictor-corrector  integration  method,  as  implemented  by  Shampine  and 
Gordon  (see  Ref.  5),  has  been  made  available  as  an  alternative  to  the  standard  second-order 
Runge-Kutta  method  used  in  BOMMP.  The  Adams  method  automatically  optimizes  the  order 
and  integration  interval  to  provide  maximum  computing  speed  for  specified  accuracy  and  round-off 
error.  This  method  provides  large  reductions  in  computing  time  when  solving  non-stiff  systems 
of  differential  equations  in  which  significant  periods  of  inactivity  occur  and  in  which  the  derivatives 
are  expensive  to  compute. 
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The  second-order  Runge-Kutta  method  may  be  used  in  all  versions.  For  the  versions  that 
include  the  Adams  integrator,  there  are  slight  differences  to  the  ‘INTEGRATION’  command. 
The  message 

ADAMS  INTEGRATOR: 

appears  first,  and  the  user  must  type  “Y”  to  specify  the  Adams  integrator,  or  any  other  character 
(including  a ‘carriage-return’)  to  specify  the  Runge-Kutta  second  order  integrator.  The  rest  of 
the  prompts  should  be  answered  as  before.  When  the  Adams  integrator  is  specified,  the  integration 
interval  is  ignored.  The  Adams  integrator  automatically  varies  its  order  and  integration  interval 
for  maximum  efficiency  within  specified  relative  and  absolute  errors.  These  may  be  specified 
by  including  them,  in  the  above  order  suitable  for  input  via  the  statement  FORMAT  (2F),  in 
a file  on  logical  unit  37  called  ADAMS. ERR.  If  this  fil^  is  not  provided,  the  values  0 0001  and 
0-0000001  respectively  are  used. 

2.1.2  DEC  FORTRAN-10  Compilation 

Minor  changes  have  been  made  to  the  FORTRAN  coding  to  allow  compilation  using  the 
F10  compiler,  which  can  produce  programs  that  run  up  to  40%  more  efficiently  in  CPU  time. 
However,  no  routine  is  presently  available  to  provide  expandable  arrays  in  F10  compiled  pro- 
grams. In  F10  versions  therefore,  there  is  a maximum  block  number  of  300  and  a maximum 
number  of  1 and  Tl,  U,  and  F blocks  of  100,  25,  and  25  respectively.  These  limits  may  be  altered 
by  changing  the  appropriate  DIMENSION  statements  in  the  files  MAIN2.FOR  or  MAIN4.FOR 
(see  Section  2.2). 


2.1.3  The  ‘RESUME’  Command 

A new  command,  ‘RESUME’,  has  been  added.  It  enables  the  continuation  of  a run  that 
has  been  interrupted  by  either  typing  an  “t”,  or  the  execution  of  a Q block. 

2.1.4  The  ‘RETAIN’  Command 

The  ‘RETAIN’  command  previously  reset  the  initial  condition  of  I,  Tl,  U,  and  Z blocks. 
It  now  also  resets  the  initial  condition  of  UO  blocks. 


2.2  Source  Files 

The  source  files  for  BOMMP  may  be  found  on  DECtape  068  by  ARL  users.  Copies  will  be 
made  available  to  others  on  request  when  a 7-track  magnetic  tape  is  supplied.  Unless  otherwise 
requested,  the  code  will  be  recorded  at  556  b.p.i.  in  6-bit  BCD  code  with  one  84-byte  record 
per  block. 

Files  with  extensions  ‘FOR’  and  ‘MAC’  are  FORTRAN  and  MACRO- 10  coded  respectively. 
Each  version  requires  the  files  CSMPA.FOR,  CSMPB.FOR,  CSMPC.FOR,  FCHECK.MAC, 
CPU.MAC,  TTYCHK.MAC,  and  DUSER.FOR.  Additional  files  required  for  individual 
versions  are  given  in  Table  2. 


TABLE  2 

For  each  version  of  BOMMP,  files  required  in  addition  to 
CSMPA.FOR,  CSMPB.FOR,  CSMPC.FOR,  FCHECK.MAC, 
CPU.MAC,  TTYCHK.MAC,  and  DUSER.FOR 


Version 

Additional  Files 

1 

MAINI.FOR, 

EXPAND. MAC 

2 

MAIN2.FOR 

3 

MAIN3.FOR, 

EXPAND. MAC,  INTEG.FOR 

4 

MAIN4.FOR, 

INTEG.FOR 

2 


2.3  Running  BOMMP  on  the  ARL  PDP-10  Computer 


2.3.1  Relocatable  Binary  Files 


For  users  on  the  ARL  PDP-10  computer,  the  relocatable  binary  files  may  be  found  on 
DECtape  137.  The  source  files,  the  compiler  used,  and  the  subsequent  ‘REL’  files  are  given 
in  Table  3.  Where  an  extension  is  not  shown,  it  is  assumed  to  be  ‘REL’. 


TABLE  3 

Relocatable  binary  files  and  their  sources 


FORTRAN  Source 

Relocatable  binary  files 

F40  Compiler  F10  Compiler 

MAIN1.FOR 

MAIN1 

— 

MAIN2.FOR 

— 

MAIN2 

MAIN3.FOR 

MAIN3 

— 

MAIN4.FOR 

— 

MAIN4 

CSMPA.FOR,  CSMPB.FOR,  CSMPC.FOR 

SUB40 

SUB10 

1NTEG.FOR 

1NT40 

INT10 

DUSER.FOR 

DUS40 

DUS10 

The  ‘REL’  versions  of  the  MACRO  routines  are  combined  in  the  file  MACROS. REL. 
The  Version  1 and  2 ‘REL’  files  have  been  combined  as  follows  (extension  ‘REL’  is  assumed 
as  above) 

BOM40  = MAIN1  + SUB40  + MACROS 
BOMMP  = MAIN2  + SUB10  + MACROS 
and  may  be  found  on  DSKB:  [1031,  1161]. 


2.3.2  Loading  the  Progrhm 

Assuming  that  the  required  ‘REL’  files  are  on  the  user’s  own  disk  area,  then  a core  image 
of  the  program  may  be  produced  in  the  file  BOMMP.SAV  for  each  version  as  follows.  The  files 
are  loaded  by  using  the  system  program  ‘LINK’. 

a)  Version  1 
. R LINK 

*BOM40  (or  MAIN1,  SUB40,  MACROS) 

♦(the  ‘REL’  versions  of  any  user-defined  subroutines) 

♦/SEARCH  DUS40 
♦BOMMP/SA  VE/GO 

EXIT 


b)  Version  2 


. R LINK 

♦BOMMP  (or  MAIN2,  SUB10,  MACROS) 

♦(the  ‘REL’  version  of  any  user-defined  subroutines) 

‘/SEARCH  DUS  10 

♦BOMMP/SAVE/GO 

EXIT 


c)  Version  3 
. R LINK 

♦MAIN3,  SUB40,  INT40,  MACROS 

*(the  ‘REL’  versions  of  any  user-defined  subroutines) 

♦/SEARCH  DUS40 
♦BOMMP/SAVE/GO 

EXIT 

d)  Version  4 
. R LINK 

♦MAIN4,  SUB  10,  INT10,  MACROS 

♦(the  ‘REL’  versions  of  any  user-defined  subroutines) 

♦/SEARCH  DUS  10 
♦BOMMP/SAVE/GO 

EXIT 

2.4  Running  BOMMP  on  Other  Machines 

In  order  to  run  BOMMP  on  other  machines,  the  following  routines,  which  are  MACRO-10 
coded,  should  be  replaced  by  equivalent  routines.  In  some  cases,  it  is  necessary  to  perform  the 
same  operation,  while  in  others,  it  is  sufficient  to  provide  a dummy  routine  thereby  removing 
the  appropriate  facility.  The  routines  are  described  below  (see  Ref.  1 for  details  on  facilities 
provided  by  routines)  and,  where  a dummy  replacement  routine  is  possible,  the  replacement 
procedure  is  given. 

a)  SUBROUTINE  EXPAND  (NS1ZE,  ARRAY) 

Description  — Expands  the  size  of  array  ARRAY  to  NSIZE  elements. 

Replacement  — Use  only  Version  2 or  4. 
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b)  LOGICAL  FUNCTION  FCHECK  (NAME) 

Description  — Searches  the  user’s  disk  area  for  the  filename  in  NAME.  FCHECK  is 
.TRUE,  if  the  file  is  found. 

Replacement  — For  the  first  entry,  set  FCHECK  to  .FALSE,  if  running  interactively, 
or  to  .TRUE,  if  running  non-interactively  (Batch  mode).  On  sub- 
sequent entries,  set  FCHECK  to  .TRUE.  . 

c)  SUBROUTINE  CPU  (TCPU) 

Description  For  the  initial  entry,  TCPU  is  the  CPU  time  in  seconds  since  the 
beginning  of  the  job.  On  subsequent  entries,  TCPU  is  the  CPU  time 
since  the  previous  entry. 

Replacement  — Set  TCPU  equal  to  zero. 

d)  LOGICAL  FUNCTION  TTYCHK  (DUMMY) 

Description  — DUMMY  is  a dummy  parameter.  TTYCHK  is  .TRUE,  if  an  “f” 
(ASCII  code  136)  has  been  typed  at  the  terminal  to  halt  model  execution. 

Replacement  — Set  TTYCHK  to  .FALSE.  . 

e)  SUBROUTINE  CHARSH  (ANS,  X,  I) 

Description  — CHARSH  operates  on  an  ascii  string  of  five  characters  contained  in  X, 
and  rotates  or  shifts  them  depending  on  I.  The  result  is  returned  in 
ANS.  For  I positive,  rotation  to  the  left  occurs.  For  I negative,  shifting 
to  the  right  occurs  with  null  characters  entered  from  the  left. 

Replacement  — Not  possible  without  extensive  recording. 

f)  SUBROUTINE  RUN  (PRONAM) 

Description  — Exits  from  BOMMP  and  then  runs  the  program  in  PRONAM. 

Replacement  — Print  on  terminal  that  the  ‘PROGRAM’  command  is  invalid  (do  not 
use  the  ‘PROGRAM’  command).  Control  automatically  returns  to  the 
command  mode. 


When  using  the  Adams  integrator,  some  modifications  are  necessary  to  the  FORTRAN 
coded  subroutine  STEP  in  the  file  INTEG.FOR.  For  the  particular  computer  being  used,  the 
round-off  error  U should  be  determined:  it  is  defined  as  the  smallest  value  such  that  (1  -0  + U) 
is  greater  than  10.  The  variables  TWOU  (i.e.  2 x U)  and  FOURU  (i.e.  4 x U)  should  be 
set  in  the  appropriate  DATA  statement.  It  should  be  noted  that  the  subroutine  STEP  and  the 
subroutine  F,  which  is  declared  EXTERNAL  in  STEP,  have  different  parameter  lists  to  those 
in  Reference  3. 

In  BOMMP,  the  DEC  FORTRAN  statements ‘OPEN’ and  ‘CLOSE’,  which  perform  basic 
file  handling  operations,  may  need  to  be  replaced. 


3.  THE  OUPUT  PROGRAM  ‘TRANS* 

3.1  Improvements 

Most  of  the  improvements  described  below  have  been  implemented  by  using  new  commands. 
These  are  listed  in  Table  4 together  with  a reference  to  the  appropriate  section  in  this  document. 
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TABLE  4 

New  TRANS  Commands 


Command 

Reference 

BARS 

3.1.11 

CYCLE 

3.1.2 

LIST 

3.1.4 

PRKONSTANT 

3.1.5 

REPEAT 

3.1.9 

SPACE 

3.1.10 

XBLOCK 

3.1.7 

3.1.1  DEC  FORTRAN-10  Compilation 

As  with  BOMMP,  minor  FORTRAN  code  changes  have  been  made  to  allow  compilation 
using  the  F10  compiler.  Because  of  the  resulting  absence  of  the  expandable  array  facility,  when 
using  this  compiler,  there  is  a limit  of  300  blocks  that  can  be  processed  for  each  input  file.  This 
limit  can  be  altered  by  changing  the  appropriate  DIMENSION  statements  in  the  main  program 
of  the  file  TRANSA.FOR  (see  Section  3.2). 


3.1.2  Non-Interactive  Running  and  Batch  Processing 

As  with  BOMMP,  TRANS  may  now  be  run  non-interactively,  during  On-Line  or  Batch 
processing.  When  running  TRANS  interactively,  all  commands  and  answers  to  prompts  necessary 
for  similar  processing  of  other  input  files  are  written  into  the  file  TRANS. OUT.  By  renaming 
this  file  TRANS.IN,  subsequent  input  files  are  processed  non-interactively:  instructions  are 
read  from  the  file  TRANS.IN.  The  renaming  may  be  achieved  automatically  while  still  running 
TRANS  (i.e.  before  using  the  ‘EXIT’  command)  by  using  a new  command,  ‘CYCLE’.  Non- 
interactive  running  requires  only  the  filename  to  be  specified,  with  a ‘carriage-return’  terminating 
program  execution.  To  revert  to  interactive  running,  the  file  TRANS.IN  must  be  renamed  or 
deleted  from  the  user’s  disk  area.  Because  not  all  answers  typed  by  the  user  are  written  into 
TRANS.OUT  (e.g.  when  using  the  ‘SCALE’  or  ‘LABEL’  commands),  the  user  is  advised  against 
creating  his  own  TRANS.IN  file  independently  of  running  TRANS. 


3.1.3  Preservation  of  Output  Files 

Assuming  an  input  file  named  A. DAT,  the  output  files  created  by  TRANS  for  tabular  or 
graphical  output  are  A. COL,  A. PPL,  and  A.PLT.  If  an  input  file  with  the  same  name  (i.e.  A. DAT) 
were  subsequently  processed,  these  output  files  would  be  overwritten  or  deleted  when  using  the 
old  version  of  TRANS.  With  the  new  version,  the  extension  name  is  modified  to  allow  preser- 
vation of  up  to  101  files  for  each  type  of  output.  For  the  above  example,  these  files  are: 


A. COL, 

A.COn 

(for  n = 0 . . 

• 9), 

A.Cnn  (for  nn  = 10  . . 

. .99) 

A. PPL. 

A.LPn 

(for  n = 0 . . 

.9), 

A.Lnn  (for  nn  = 10 . 

. .99) 

A.PLT, 

A.PLn 

(for  n = 0 . 

.9), 

A.Pnn  (for  nn  = 10  . . 

.99) 

For  more  than  101  files,  the  output  is  appended  to  the  last  appropriate  file  (e.g.  A.C99).  Following 
the  ‘PRC’,  ‘PRP’,  ‘PLS’,  and  ‘PLO’  commands,  the  output  file  is  specified  on  the  terminal  if 
the  particular  output  filename  has  been  modified  as  above  (for  non-interactive  runs  as  well). 
For  example,  after  a ‘PLS’  or  ‘PLO’  command  when  the  files  A.PLT  . . . A PI  1 already  exist  on 
the  user’s  disk  area,  the  following  message  is  typed: 

[PLS/O  output,  for  this  run,  going  to  DSK:A  .PI 2] 

If  all  the  files  A.PLT  . . . A.P99  exist,  then  the  message  typed  is: 

[PLS/O  output,  for  this  run,  appended  to  DSK.A  .P99] 


6 


3.1.4  Specifying  Block  Numbers 

When  specifying  block  numbers  in  response  to  the  prompt: 

BLKS 

the  user  must  now  type  “AV”  (previously  “A”)  to  specify  all  blocks.  In  all  instances,  except 
when  specifying  output  for  the  line  printer  using  the  ‘PRC’  or  ‘PRP’  commands  (see  Sections 

3.1.5  and  3.1.6),  typing  “A”  now  specifies  only  those  blocks  whose  value  varies. 

The  user  may  also  specify  a list,  which  consists  of  up  to  30  block  numbers,  by  typing  “Ln”, 
where  n = 1 . . . 9.  The  ‘LIST’  command  is  used  to  define  the  composition  of  each  list.  All  block 
number  lists  are  stored  in  the  file  TRANS. BLK  and  are  therefore  available  to  further  runs 
of  TRANS. 

3.1.5  Revised  Operation  of  ‘PRCOLUMN’  Command 

When  specifying  the  blocks  to  be  printed  on  the  line  printer,  by  typing  "A”,  the  constant 
values  of  the  non-varying  blocks  are  first  printed  followed  by  the  varying  block  values  in  tabular 
form.  By  using  the  ‘PRKONSTANT’  command,  only  the  constant  values  of  the  non-varying 
blocks  are  output  for  printing  on  the  line  printer. 

3.1.6  Revised  Operation  of  ‘PRPLOT’  Command 

When  specifying  the  blocks  to  be  plotted  on  the  line  printer,  by  typing  “A”,  only  the  varying 
blocks  are  plotted.  Any  non-varying  blocks  are  printed  before  the  plots  as  if  the  ‘PRKONSTANT’ 
command  had  first  been  used. 

3.1.7  Cross  Plots 

The  ‘XBLOCK’  command  allows  the  X axis  block  to  be  redefined  for  incremental  plots,  thus 
enabling  ‘cross  plots'  of  dependant  variables. 

3.1.8  Revised  Form  of  ‘Overlay’  Plots 

For  ‘overlay’  plots,  curves  previously  labelled  by  different  symbols  are  now  identified  by 
different  broken  lines.  This  was  found  to  give  a less  cluttered  appearance  to  the  plots. 

3.1.9  Plotting  Data  from  Different  Input  Files  on  the  Same  Graphs 

The  ‘REPEAT’  command  is  equivalent  to  the  ‘CYCLE’  command  except  that,  for  ‘strip’ 
plots,  subsequently  processed  files  are  plotted  on  the  same  graph(s).  The  different  curves  are 
identified  by  different  broken  lines.  In  the  key  for  each  graph,  the  broken  lines  are  labelled  with 
the  titles  (as  they  appear  in  the  input  files)  and  the  input  filenames.  The  title  at  the  base  of  each 
graph  is  that  stored  by  TRANS  on  executing  the  first  ‘GOE’  command.  The  user  is  advised 
against  using  ‘overlay’  plots  when  using  the  ‘REPEAT’  command. 

3.1.10  Spacing  between  Incremental  Plots 

Incremental  plots  are  suitably  spaced  on  8-5  inch  folding  paper  for  the  pen  initially  set  at 
a fold.  When  plotting  on  non-folding  paper,  unnecessary  spacing  can  be  removed  by  using  the 
‘SPACE’  command. 

3.1.11  Bars  on  Incremental  Plots 

The  ‘BARS’  command  results  in  the  placement  of  bars  on  the  incremental  plot  curves  at 
specified  equal  independent  variable  intervals.  The  bars  are  short  lines  perpendicular  to  the  curve 
and  they  may  be  used  to  provide  a third  demension,  e.g.  on  a ‘cross  plot’  of  X and  Y coordinates 
of  an  aircraft,  the  bar  spacing  will  indicate  speed  if  the  independant  variable  is  time. 

3.2  Source  Files 

The  source  files  for  TRANS  may  be  found  on  DECtape  068.  Files  with  extensions  ‘FOR’ 


7 


and  ‘MAC’  are  FORTRAN  and  MACRO  coded  respectively.  The  files  are  TRANSA.FOR, 
TRANSB.FOR,  FCHECK.MAC,  and  SYMBOL.MAC. 

3.3  Running  TRANS  on  the  ARL  PDP-10  Computer 

The  source  files  are  first  compiled,  using  either  the  F10  or  F40  compiler  for  the  FORTRAN 
files.  Then,  assuming  that  the  required  ‘REL’  files  are  on  the  user’s  own  disk  area,  a core  image 
of  the  program  may  be  produced  in  the  file  TRANS. SAV  as  follows.  The  files  are  loaded  using 
the  system  program  ‘LINK’. 

. R LINK 

*TRANSA,  TRANSB,  FCHECK,  SYMBOL 

♦TRANS/SAVE/GO 

EXIT 


Generally,  the  above  procedure  is  unnecessary,  since  a copy  of  TRANS. SAV  may  be  found 
on  DSKB:[103I,  1161]. 

3.4  Running  TRANS  on  Other  Machines 

In  order  to  run  TRANS  on  other  machines,  the  machine  coded  routines  FCHECK, 
SYMBOL,  and  PLOT  should  be  replaced  by  equivalent  routines.  A description  and  suggested 
replacement  of  FCHECK  is  given  in  Section  2.4.  The  routines  SYMBOL  and  PLOT  are  used 
extensively  in  generating  incremental  plotter  output:  replacement  routines  are  therefore  essential 
if  incremental  plots  are  required.  The  descriptions  of  these  two  routines  are  given  below. 

a)  SUBROUTINE  SYMBOL  (LU,  X,  Y,  HEIGHT,  ASCII,  THETA,  N) 

The  subroutine  SYMBOL  plots  the  ascii  string  of  N characters  contained  in  ASCII  at 
the  point  X,  Y.  LU  is  the  logical  unit  on  which  the  plot  file  is  being  written.  HEIGHT  is  the 
height  (in  inches)  of  the  characters  and  THETA  is  the  angle  (in  degrees)  to  the  x-axis  at  which 
the  string  is  to  be  plotted. 

b)  SUBROUTINE  PLOT  (LU,  X,  Y,  N) 

The  subroutine  PLOT  controls  the  pen  on  the  incremental  plotter.  LU  is  the  same  as  in 
SYMBOL.  X and  Y are  coordinates  if  they  are  real  variables,  or  the  number  of  0-01  inch  steps 
moved  in  the  X and  Y directions  if  they  are  integer.  X and  Y must  be  real  for  N = 1 or  2. 
For  N = 1,  X and  Y are  scaling  factors.  For  N = 2,  a new  origin  is  defined  such  that  X and  Y 
are  the  coordinates  of  the  current  position  of  the  pen  with  respect  to  the  new  origin.  For  N = 3, 
the  pen  is  lifted,  moved  to  X,  Y and  dropped.  For  N = 4,  the  pen  is  dropped  and  moved  to  X,  Y. 
For  N = 5,  the  pen  is  lifted,  moved  to  X,  Y but  not  dropped. 

As  with  BOMMP,  the  DEC  FORTRAN  statements  ‘OPEN’  and  ‘CLOSE’  may  need  to  be 
replaced. 

4.  CONCLUDING  REMARKS 

This  postscript  to  previous  publications  on  the  simulation  language  CSMP-10(ARL) 
consists  of  a description  of  improvements,  details  of  the  location  of  source  files  and  their  use 
on  the  ARL  computer,  and  suggestions  on  how  to  implement  the  language  on  other  computers. 
The  improvements  significantly  increase  the  efficiency  of  the  language  in  handling  simulation 
problems.  DEC  FORTRAN-10  compilation  and  an  alternative  integration  method  that  is  very 
efficient  for  certain  systems  make  possible  large  savings  in  CPU  time.  The  many  improvements 
to  the  output  program  enable  it  to  be  used  with  much  greater  flexibility,  especially  when  handling 
several  output  files  concurrently,  e.g.  plotting  data  from  different  files  on  the  same  graph.  The 
aircraft  arresting  gear  problem  of  Reference  4 is  used  in  a comprehensive  demonstration  of  the 
new  operating  procedures  and  improvements  to  the  simulation  language. 
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APPENDIX  A 


Use  of  CSMP-10(ARL)  for  the  Aircraft  Arresting  Gear  Problem 

A.1  Description  of  Problem 

(1)  Physical  diagram 


Cable  2 

spring  constant 


Cable  1 

spring  constant 


(2)  Equations 

m$  y 3 = /* 2 — ft> 

MO)  = MO)  = 0 

ms  ys  = 2/*i  —fa 

MO)  = MO)  = o 

mix  = — 2/*i  sin  8 

x (0)  = 0,  x(0)  is  variable 

fkS  = k 2 (^2  — J's) 

if  ys  < ys 

= 0 

if  I’s  > ys 

fki  = ki  (yi  — 2ys) 

if  yi  > 2ys 

= 0 

if  yi  2ys 

fD  = f(ya)  ys2 

where  f(ys)  is  defined  in  (4) 

yi  = (x2  + h2)*  - h 

sin  8 = x/(h  + >>i) 

(3)  Constants 

mi  = 1400  slug 
m 2 = 45  -28  slug 
ms  = 20  slug 
*i  = 4550  lb/ft 
k2  - 25300  lb/ft 
h = 125  ft 


(4)  Function  f (>>3) 

Defined  by  linear  interpolation  of  the 
following  coordinates : 


ya 

/(ya) 

(ft) 

(Ib/(ft/s)*) 

0 

8-33 

30 

4-0 

60 

1-6 

120 

5-2 

150 

5-2 

180 

6-6 

210 

8-3 

240 

10-7 

270 

16-0 

282 

21  -0 

294 

28-0 

306 

410 

312 

50-0 

324 

90-0 

999 

90-0 

(5)  Output 

a)  At  x(0)  = 290  ft/s  using  both  Adams  and  Runge-Kutta  integrators  for  time 
t = 0 to  10s, 

(i)  tabular  output  of  £3,  ^3, /n,  y2,  x,  x in  Is  intervals; 

(ii)  graphical  output*  of  x,  x,  y3  vs.  t;  and 

(iii)  graphical  output*  of  x vs.  x. 

b)  Graphical  output  of  x,  x,  y3  vs.  t for  t = 0 to  10s  for  x0  = 290,  250,  210  ft/s 
(on  same  graph)  using  Adams  integrator. 


* Since  graphs  are  almost  identical  for  runs  with  either  integrator  used  (assuming  the  same 
solution  accuracy),  both  graphs  are  not  shown  below. 


(6)  Block  diagram  (unless  shown,  integrator  initial  conditions  are  zero;  see  Fig.  1 of  Ref.  1 
for  notation) 


1.  Block  10  output  = ya  — yz  if  < yz 

= 0 if  )3  3*  yz 

2.  Block  26  output  = yi  — lyz  if  yi  > 2yz 

= 0 if  >i  < 2yz 


A.2  Model  Input  File 


CONFIGURATIONS  PARAMETERS 


X 

DOT(0) 

* 291 

12  2.ffl«E+f 1 

2 

X 

7 

7 

>3  -2.5300E+04 

3 

F 

8 

14  4.5280E+0I 

4 

X 

2 

3 

; fi 

19  1 .259IE+I2 

S 

♦ 

11 

-4 

28  9.10I0E+03 

6 

/ 

5 

12 

;T3  iot  dot 

31  1 .4M0E+03 

7 

I 

6 

; Y3  DOT 

32  2.900IE+I2 

8 

I 

7 

; y3 

9 

♦ 

8 

-18 

FUNCTIONS 

II 

P 

9 

3 

11 

X 

11 

13 

; FK2 

f .I0I0E-0I 

I.3300E+I0 

12 

K 

; M3 

3.lfffEtf 1 

4.0I00E+0I 

13 

K 

; -K2 

6.llffE+f 1 

1.6I00E+0I 

14 

K 

; M2 

1 .2fffE+f2 

5.2I00E+0I 

15 

♦ 

-11 

27 

1 . 5###E + #2 

5.200IE+I0 

16 

/ 

15 

14 

;Y2  dot  dot 

1.8IME+I2 

6.600IE+I0 

17 

I 

16 

; T2  IOT 

2.10IIE+I2 

8.3I0IE+I0 

18 

I 

17 

; Y2 

2.4II0E+I2 

1 .I70IE+I1 

19 

K 

; h 

2.70ffE+02 

1 .600IE+1 1 

23 

U1 

33 

19 

; Y1 

2.82ME+I2 

2.100IE+I1 

24 

UO 

23 

;SIN( THETA) 

2.94IIE+02 

2.800IE+I1 

25 

♦ 

-18 

-18 

23 

3.I69IE+I2 

4. 100IE+I 1 

26 

N 

25 

3.12ME+I2 

5.I00IE+I1 

27 

X 

26 

28 

; 2FKI 

3.2400E+02 

9.000IE+01 

28 

K 

; 2Ki 

9.9900E+I2 

9.I0I0E+01 

29 

- 

35 

31 

K 

; «i 

31 

/ 

29 

31 

;x  DOT  DOT 

32 

I 

31 

; x IOT 

33 

I 

32 

; x 

35 

X 

27 

24 

;2FKt*SIN 

L_ 


I 


oooooooo 


A3  User-defined  Subroutine 


SUBROUTINE  USER1 (II,C,NTIX,PAI) 

USER-DEFINED  SUBROUTINE  FOR 
THE  ARRESTIN6  BEAR  PROBLEM 

IT  CALCULATES 

T1«SOIT<X*X*H*H)-H 
AS  ITS  PRINCIPLE  OUTPUT,  AN I 
SIN(THETA)*X/(I+Y1 ) 

AS  AN  ADDITIONAL  OUTPUT  VIA  A UO  BLOCK. 

DIMENSION  C(1),NTRX(1),PAR(1) 

COMMON  /TEST/TEST (») 

INTE6ER  TEST 
DATA  IIUOt/24/ 

C FOR  INITIAL  EXECUTION  TEST(5)»I 

IF ( TEST < 5 > . NE . 1 ) SO  TO  1* 

C 6ET  INPUT  BLOCK  NUMBERS 

IB1=NTRX(5*(IB-1 )*2) 
ID2»MTRX(5*(IB-1)+3) 

C 6ET  INPUT  VALUES 

11  X=C(IB1 ) 

H*C(IB2) 

C CALCULATE 

Y1*SORT (X*X+H*H)-H 
STHETA«X/(H+Y1) 

C SET  OUTPUT  VALUES 

C(IB)*Y1 
C< IBUOt )»STHETA 

RETURN 

END 


A.4  Creating  and  Running  Modelling  Program 

(1)  Creating  stored  program  using  Version  4 of  BOM MP 


•C0HPILE/F1 f/NEU  EXIISER.FOR 

FORTRAN i EXUSER 

USER1 

•R  LINK 

♦MAINA ,SUB1 f , INTI » f MACROS 
♦EXUSER 

♦/SEARCH  DUSt# 
♦ARSTR/SAVE/60 

EXIT 


(2)  Running  program  with  Adams  integrator  for  x0  = 290,  250,  210 
(uses  default  Adams  integration  parameters) 


.RB  AftSTR 

MAI  BLK  NO.  * 3M 

MAX  NO.  OF i I S T1  BLXS,  U DLKS,  F ILKS  « If I, 23, 25 
INTEGRATORS: 

RUN6E-KUTTA  2ND  ORDER  C DEFAULT] 

ADANS  VARIABLE  STEF  (M.K. GORDON  IMPLEMENTATION) 

♦1068: ARSTRtI 

MODEL  I/P  FRON  LOGS: ARSTR. NOD 
♦LOG?  i OUT  UO 

BLOCK  0/P  TO  LOI9IOUTI  .DAT 


♦COM 


CONFIGURATIONS  i 
X DOT(f)  « 2P0 

ILK  TYPE  II  B2  B3 

2 X 7 7 

3 F 8 

4 X tO 

♦PAR 

PARAMETERS  t 

ILK  PI  P2  P3 

1 12  2.0000E+01 

13  -2.330IE+04tO 

♦FUN 

FUNCTIONS  i 
BLK  NO.  3 
COORD  PAIRS  : 


I .I090E+0I  8.3J00E+0I 

3.ffffE+f 1 4.0I00E+0HO 


RUN  CPU  TINE  : 15.78  SEC 
♦TIT 


TITLE  (LIMIT  if  CHRS) 

X DOT (I)  « 251 

♦PAR 

PARAMETERS  i 

BLK,  PI,  P2,  P3 
32,23# 

♦LOG9:OUT2«-0 

BLOCK  O/P  TO  LOS9:OUT2  .DAT 
♦DOE 

♦♦  RUNNIN6  ♦• 

RUN  CPU  TIME  t 15.07  SEC. 

♦TIT 

TITLE  (LIMIT  if  CHRS) 

X DOT(f)  « 211 

♦PAR 

PARAMETERS  t 

BLK,  PI,  P2,  PJ 
32,210 

♦LQG9sOUT3+l 

BLOCK  O/P  TO  LOS9IOUT3  .RAT 
♦60E 

♦•  RUNNING  ♦• 

RUN  CPU  TIME  i 13.64  SEC. 

♦EXI 

END  OF  EXECUTION 

CPU  TINEi  49.18  ELAPSED  TIME:  7i2t.6U 
El 


I 


(3)  Preparing  for  further  runs  that  are  partly  non-interactive 


.1  PIP 

*IOHMP.INb)OHNP.0UT/T 

*»c 

.TYPE  B0HHP.1N 
L068:  ARSTR«-I 
L069:OUTHO 
CON 

PAR 

FUN 

INT 

0,10,9 

Y 

OUT 

A 

0.1001,0.1 

OOE 

TIT 

X IOT(0)  - 25f 
PAR 

32,231 

L0G9:0UT2H) 

GOE 

TIT 

X DOT(ff)  • 21# 

PAR 

32,211 

L069iOUT3+O 

GOE 

EXI 

.TECO  BOHHP.IN 

C2K  CORE! 

*SINT«fLT«f 

INT 

*3#KTtf 

♦IRAN 

»*3T« 

FUN 

HAN 

+EXW 


(4)  Running  program  with  Runge  Kutta  integrator  for  Xo  = 290 

(program  is  run  partly  non-interactive;  shows  use  of  (i)  ‘RESUME’  following  “f”  and 
(ii)  ‘RETAIN’ 

.RU  ARSTR 

MAX  BLK  NO.  = 3M 

MAX  NO.  OF:  I S T1  BLKS,  U ILKS,  F ILKS  * 1#», 25,25 
INTEGRATORS: 

RUN6E-KUTTA  2ND  ORDER  [DEFAULT] 

ADAMS  VARIADLE  STEP  (M.K. GORDON  IMPLEMENTATION) 

X DOTH ) » 29f 
*L069:0UTRKt0 

BLOCK  0/P  TO  L0G9:0UTRK.DAT 
♦INT 

INTEGN  PARAHS;  LOUER,  UPPER,  INTERVAL  =»  0,t«,f.#D5 

ADAMS  INTEGRATOR  : N 
♦OUT 

0/P  DLKS 
A 

0/P  PARAMS;  Z CHANGE  REQRD,  INTERVAL  = 


♦60E 


♦*  RUNNING  M‘ 


RUN  CPU  TINE  : 6.48  SEC. 

RUN  TERMINATED  BY  A 

♦LOO 
BLK  = 1 

O/P  = 2.6IME+0I 
BLK  - 32 

O/P  = 1 .8303E+I2 

BLK  = 

♦RES 

♦♦  RUNNING  •♦ 

RUN  CPU  TIME  : 25.94  SEC. 

♦LOO 
BLK  « 1 

O/P  > 1 . 00N0E+0I 
BLK  s 32 

O/P  * 5.9902E+I0 
BLK  > 

♦RET 

♦LOGIN  :FINSH«-M 

MODEL  O/P  TO  LOGIIsFINSH.MOD 
♦STO 

CON,  PAR,  FUN,  OR  ALL  :ALL 
♦TTTtfN 

MODEL  O/P  TO  TTY: 

♦STO 


CON,  PA*,  FUN,  OR  ALL  -.PAR 


7 

2.9794E+00 

8 

3.8383E+02 

12 

2.000IE+I 1 

13 

-2.5300E+04 

14 

4.5280E+I1 

17 

2.9736E+00 

18 

3.8386E+02 

19 

1.2500E*02 

24 

9.9015E-01 

28 

9.1O0IE+I3 

3» 

1 .4IO0E+03 

32 

5.9902E+I0 

33 

8.84O1E+02 

*EXI 

ENI  OF  EXECUTION 

CPU  TINE!  30.14  ELAPSED  TINE:  5:17.78 
EXIT 


— 1 

A.S  Running  Output  Program 

(1)  Output  of  corresponding  runs  for  xo  = 290  using  Adams  and  Runge-Kutta  integrators 
(shows  use  of ‘LIST’,  ‘PRKONSTANT’,  ‘SPACE’,  ‘XBLOCK’,  ‘BARS’,  and  ‘CYCLE’; 
files  TRANS. IN  and  TRANS. BLK,  as  well  as  tabular  and  incremental  plotter  output, 
are  presented) 

.R  PIP 

•/X+E1031 ,1161 3TRANS.SAU 
•♦C 

.RU  TRANS 

CTRANS  VERSION  DATE  24-FEB-783 
I/P  FILENAME  * 0UT1 
X DOT (•)  * 290 

I/P  FILE  RECORDED  ON  27-FEB-78  AT  11 :1t 
INTE6N  INT  = RUN  CPU  TINE  * 15.78  SEC. 

TINE  FROM  • .IMOE+IO  TO  1 .MfOE+01  IN  STEPS  OF  I .ffME-il 
•LIS 

BLOCK  LIST  NUMBER  = 1 
LIST  OF  BLOCK  NUMBERS  > 7,8,4,11,32,33 
BLOCK  LIST  NUMBER  * 2 
LIST  OF  BLOCK  NUMBERS  = 32,33,7 
BLOCK  LIST  NUMBER  » 

•TIM 

TIME  PARAHS;  LONER,  UPPER,  INTERVAL  = 1,11,1 
•PRK 
•PRC 

PRINTIN6  IN  COLUMNS  : 

BLKS 
LI 

IS  0/P  TO  TTY  REORD  : N 
•SPA 

IS  SPACING  BETUEEN  PLOTS  REORD  : N 

J 


i 


*PLS 


STRIP  PLOTS  : 

BLKS 

L2 

LENGTH  OF  AXES  IN  INCHES;  X,  Y » 3,2  « 
*60E 

*»  RUNNING  « 

*XIL 

NEW  X AXIS  ILK  NO.  = 32 
*»AR 

INDEP  VARIABLE  INTERVAL  BETWEEN  BARS  ■ I 
*PLS 

STRIP  PLOTS  s 

BLKS 

33 

LEN6TH  OF  AXES  IN  INCHES;  X,  T *5,5 
*60E 

**  RUNNING  ** 

*CTC 

I/P  FILENAME  * QUIRK 

**  RUNNING  ** 

**  RUNNING  •* 

I/P  FILENAME  = 


END  OF  EXECUTION 

CPU  TIHEt  1 A .34  ELAPSED  TINEl  3i34.A2 
EXIT 


.TYPE  TRANS. IN 

US 

TIN 

MHIE+II  1 .N00IE+I1  1.0N0NE+I0 

PRK 

PRC 

LI 

N 

SPA 

N 

PLS 

L2 


3.0  2.0 


GOE 

XIL 

32 

BAR 

1 .00I0E+00 

PLS 

33 


5.0  5.0 
GOE 
CYC 

.DELETE  TRANS. IN 
FILES  DELETED: 

TRANS. IN 
01  BLOCKS  FREED 

.TYPE  TRANS. BLK 
LI 

7 8 4 10  32  33 


L2 

32  33  7fl#f|00f00«0|0» 

»IINI»0lffl« 

.R  PRINT 


11:34*0UT1. COL/P. OUTRK. 

CQL/PN 

FILES  PRINTED 

RED. 

PANES 

DSKtOUTI .COL 

Rf  1 

2 

DSK:OUTRK.COL 

EXIT 

R02 

2 

.R  PLOTQ 

1 1 :37*0UTRK.PLT$ 

FILES  PLOTTED 

RES. 

FEET 

LIMITS 

OUTRK  . PLT 
EXIT 

RN1 

2 

X(3.28,F1/7M) ,T(-5.32"f5"> 

Tabular  output  file  OUT1.COL 

(uses  Adams  integrator) 


X DOT (•)  • 29# 


RECORDED  ON  27-Feb-78  AT  11 sli  INTEBN  INT  » «.8##fE+l#  RUN  CPU  TIUE  « 15.78  SEC. 
NON-VARYING  BLOCK  VALUES 


ILK  NUNBER  = 12 

13 

14 

19 

21 

3# 

N3 

-K2 

N2 

H 

2K1 

HI 

2.#fffE+#1 

-2.3388E+84 

4.5288E+81 

1.25»#E+#2 

9.1#0»E+»3 

1 .4MIE+I3 

X DOTH) 

RECORDED  ON 

= 29» 

27-Feb-78  AT  11:14  INTESN  INT  * 8 

.0I00E+0I  RUN  CPU  TINE 

= 15.78  SEC, 

ILK  NUNBER  » 

' 7 

8 

4 

18 

32 

33 

TINE 

Y3  DOT 

Y3 

FD 

Y2 

X DOT 

X 

l.flHEtlf 

I .9f#fE+ff 

I . VffVE+ff 

f .IfffE+ff 

f .f«8fE+0f 

2.9»##E+*2 

0.I00IE+I0 

1 JIHEtfl 

1 .8882E+f2 

8.7538E+81 

3.7942E+I4 

8.911 4E 1 

2.79221+92 

2.8151E+02 

2.MIIEt|| 

1.8512E+82 

2.U153E+82 

8.6415E+84 

2.I493E+I2 

2.22I5E402 

5.296AE+02 

3.ffffE+fl 

6.9323E+8 1 

2.9878E+I2 

1 .2532E+I5 

2.9539E+02 

1.5190E+02 

7.1871E*#2 

4.9ffiE+ff 

3.8738E+81 

3.4I32E+82 

1.35ME+85 

3.4559E+02 

5.I044E+81 

8.21 15E+02 

5.f#«fE*ff 

7. 4788E+88 

3.A355E+82 

5.3848E*«3 

3.4375E+U2 

1 . 1254E+I1 

8.4378E+82 

4.ffffE+ff 

4.7752E+08 

3.6898E+I2 

2.I522E+I3 

3.6898E+82 

9.5154E+I0 

8.5409E+82 

7.ff#fEt«# 

4. 1 5I9E+M 

3.7334E+82 

1.55I7E+I3 

3.7348E+92 

8.2895E+I0 

8.4294E+82 

8.MffE+if 

3.6A11E+M 

3.7723E+I2 

1.2843E+83 

3. 7728E+82 

7.3470E+08 

8.7I7AE+02 

9.#mE«ll 

3.2823E+II 

3.8878E+I2 

9.496'3E*I2 

3.8874E+82 

6.59B8E+00 

8.7772E+02 

1 .ffffE+f 1 

2.9648E+M 

3.8383E+82 

7.91 11E+82 

3.8384E+I2 

5.9899E*## 

8.84»#E*#2 

Tabular  output  file  OUTRK.COL 

(uses  Runge-Kutta  integrator) 


X DOT  (f ) - 291 


RECORDED  ON  27-F»b-7l  AT  11:25  IRTE6N  INT  ■ 5.0I00E-03  RUN  CPU  TIRE  * 25.94  SEC. 
R0N-VARY1N6  BLOCK  VALUES 


ILK  NUNDER 


12 

N3 


13 

-K2 


14 

N2 


19 

H 


28 

2K1 


3f 

HI 


2.HME+I1  -2.53ME+I4  4.528*E*I1  1.25I0E+02  9.1I00E+03  1.400IE+I3 


X DOT ( « ) = 29* 


RECORDED  ON 

27-F*b-78  AT  11:25  INTE6N  INT  » 5.M*«E-f3  RUN  CPU  TIME 

= 25.94  SEC, 

ILK  NUMBER  > 

» 7 

8 

4 

18 

32 

33 

TINE 

Y3  DOT 

Y3 

FD 

Y2 

X DOT 

X 

I .ffffE+ff 

f .ffffE+ff 

«.ff«*E+** 

I.HffEtll 

I .M0IE+I0 

2.900*E*02 

0.0090E+00 

l.fflfE+fl 

1.*82*E+*2 

8.7526E+1 1 

3.8*64E+*4 

8.91 1 3E+I1 

2.7022E+I2 

2.8151E+02 

2.f«f«E+lf 

1 .*5«8E  + f2 

2.1 154E+I2 

8.6352E+I4 

2.I493E+I2 

2.2205E+I2 

5.2964E+I2 

3.IIIIEtfl 

6.9I33E+I1 

2.9*7fE+#2 

1 . 24  27E  +*5 

2.9559E+I2 

1 .5190E+02 

7.187IE+I2 

4.f*f*E«ff 

3.8362E+I1 

3.4I32E*#2 

1 .3383E+I5 

3.4559E+B2 

5.0047EHM 

8.2115E+02 

5.lf#fE+ff 

7.763IE+II 

3.6355E+I2 

5.4244E+I3 

3.6375E+02 

1 . 1255E+01 

3.4378E+02 

i.llflEtll 

4.7654E+II 

3.689IE+I2 

2.B44fE+*3 

3.6898E+I2 

9.5161E+00 

8.54I9E+02 

7.Mf*E+ff 

4.1396E+M 

3.7334E+I2 

1 . 5423E+*3 

3.734IE+I2 

8.290IE+I0 

8.62?4E*02 

8.ffffE+0f 

3.6623E+II 

3.7723E+*2 

1 . 2871 E+#3 

3.7728E+I2 

7.3473E+0I 

8.7074E+I2 

9.ffffE*if 

3.2852E+M 

3.8*7*E+*2 

9.71 33E+I2 

3.8#74E+»2 

4.5991 E+0* 

B.7772E+I2 

1 .IIIIE+11 

2.9794E+** 

3.8383E+*2 

7.9893E+I2 

3.B384E+I2 

5.99I2E+00 

8.8401E+02 

Li  1 


A A 


X DOT (0)  = 290 


. 


(2)  Plotting  three  runs  performed  in  Section  (2)  of  Appendix  A A on  same  graph 
(shows  use  of  ‘REPEAT’  and  preservation  of  output  files) 


.RU  TRANS 

i 

CTRANS  VERSION  DATE  24-FEI-78D 
I/P  FILENANE  = 0UT1 
X DOT (•)  * 2*0 

I/P  FILE  RECORDED  ON  27-FEB-78  AT  11(14 

INTE6N  INT  - f.OMfE+M;  RUN  CPU  TINE  = IS.78  SEC. 

TINE  FROM  O.OeiOE+fO  TO  1.ff«fE*01  IN  STEPS  OF  l.fl0fE-0l 
♦PIS 

CPLS/0  OUTPUT,  FOR  THIS  RUN,  GOING  TO  DSKsOUTI  .PLU 

STRIP  PLOTS  s 

BLKS 

L2 

LEN6TH  OF  AXES  IN  INCHES;  X,  Y = 5,2 
♦LAB 

BLK  NO.  -1  DENOTES  INIEP  VARIABLE 

IS  LABELLIN6  TO  BE  REAS  FROM  DSX  s N 

IS  TTY  LISTING  REQRD  : N 

ARE  HODIFICATIONS  REQRD  : Y 

BLK  <•  FOR  TITLE),  LINE,  TEXT 

0,1,  AIRCRAFT  ARREST  INS  GEAR  PROBLEM 

*8DE 

*•  RUNNING  ** 


♦REP 

I/P  FILENAME  * OUT2 


♦•  RUNNING  ♦♦ 


I/P  FILENAME  * 0UT3 


♦*  RUNNING  ** 


I/P  FILENAME  = 


END  OF  EXECUTION 

CPU  TINEs  8.96  ELAPSED  TINEi  3i9.49 
EXIT 

.TYPE  TRANS. IN 

PLS 

L2 

S.9  2.9 
LAB 
GOE 
REP 

.R  PLOTQ 


11:41* OUT I.PLIt 

FILES  PLOTTEI  REI.  FEET  LINITS 

OUT1  .PLN  Rfl  2 X(l.5N'’f  1 '9.24*),  Y(-5.  J2".5" ) 

r w * r 


Graphical  output  file  OUT1.PLO 

(a)  ‘REPEAT’  key  (appears  to  right  of  strip  plots  on  incremental  plotter  output) 


REPeot 

l 

5 

© 

x; 

1 

— i 0UT1 

X 

DOT (0)  = 290 

»_ 

-h  0UT2 

X 

DOT (0)  = 250 

i 

— i 0UT3 

X 

DOT (0)  = 210 
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